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1. Introduction 



The correspondence between the critical phenomena on flat and random lattices has 
been very useful in various problems related to random geometry. Especially interesting 
from this point of view are the 0{n) model and the q'-state Potts model 
which allow a neat geometrical interpretation in terms of a gas of loops or clusters with 
fugacity depending on the continuous parameter n or q. Considered on a flat lattice, these 
statistical models can be mapped onto the 6- vertex model whose infrared behavior is 
that of a compactified boson. The continuous parameters n and q are then related to the 
compactification radius. The mapping , known as Coulomb gas picture, also prescribes 
to install a background electric charge "at infinity" . 

On an irregular lattice, the background electric charge, which is in fact coupled to the 
local curvature, should be distributed all over the lattice, which can be achieved only by 
introducing nonlocal operators at all points with curvature. The mapping does not hold 
any more, which means that, when considered on a random lattice, the 6- vertex model 
describes a new type of critical phenomena, in which the local curvature of the lattice is 
involved. 

More concretely, let us consider the representation fl^ of the 6-vertex model as a gas 
of oriented loops. The fugacity of each loop is given by a phase factor, the phase being 
proportional to the geodesic curvature along the loop. On the flat lattice, the geodesic 
curvature is always ±2n and all loops have the same fugacity. On the contrary, on a 
lattice with defects, the geodesic curvature can take any value between n/2 and oo, and 
the fugacity of the loop depends on the local curvature of the lattice. In this way, the 
6-vertex model on a random lattice allows to study extended objects coupled to the local 
curvature, unlike all previously considered statistical models. 

The 6-vertex model on a random lattice is also interesting from the string theory point 
of view, because it gives a dual definition of the compactified c = 1 string, in which the 
vortices can be constructed explicitly as local operators. As such, it was considered first by 
P. Ginsparg [T2\, who formulated the corresponding matrix model. The relation between 
the 6-vertex model on a random lattice and the c = 1 string theory has been further 
elucidated by S. Dalley |T^, who in particular identified the compactification radius as a 
function of the vertex coupling. 

In this paper we present the exact solution of the 6v model defined on a planar 
random lattice of any finite size. The solution is obtained using the equivalence with a 
special hermitian matrix model, which is a generalization of the 0{2) matrix model 
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We will calculate the disk amplitude and investigate its scaling behavior. We will evaluate 
the fractal dimension of the boundary of the random surface as well as the average number 
of loops per unit area. We will also calculate the dimensions of the m- vortex operators 
(m > 1), which in the loop gas representation are the sources of m equally oriented vortex 
lines. Our results agree with the correspondence between the 6- vertex model on a random 
surface and the compactified c — 1 string theory, where vortices represent discontinuities 
on the world sheet. 

The paper is organized as follows. In Sect. 2 we remind the definition of the 6v 
model on a flat and random lattices as well as the correspondence with a compactified 
boson, respectively with the compactified c = 1 string theory. The model defined on a 
random lattice depends on two coupling constants: the vertex coupling A G [0, 1] and 
the " cosmological constant" 1/T coupled to the number of squares of the lattice. Then 
we formulate the partition function of the model as a matrix integral, making one step 



further than the original definition by P. Ginsparg [12|. In Sect. 3 we formulate and solve 
the saddle point equations for the resolvent of the random matrix for all allowed values 
of the two couplings. In Sect. 4 we analyze the the scaling limit of the solution near 
the critical coupling where the size of the random lattice explodes. In Sect. 5 we give 
a discussion and a summary of the results. The three appendices contain the proof of 
the exact correspondence between the 6v matrix model and the compactified c = 1 string 
theory (A), the technical details of the calculation (B), and the solution of the 0{2) matrix 
model obtained as the A ^ limit of our general solution (C). 

While this manuscript was being prepared, we received the interesting paper by Paul 
Zinn- Justin |T^, in which some of our results concerning the critical behavior are derived 
using a different technique. 



2. The 6-vertex model on a dynamical random lattice 

2.1. Review of the 6-vertex model on a flat square lattice: F -model, loop gas, vortices, 
mapping to a compactified boson. 



The degrees of freedom of the 6-vertex model on a square lattice represent arrows 
associated with the links of the lattice. The admitted configurations of arrows do not 
contain sources and sinks, so that there is an equal number of arrows that begin and end 
at each vertex. The partition function is defined as the sum over all arrow configurations, 
with Boltzmann factor, which is a product of the statistical weights associated with the 
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vertices of the lattice. On the dual lattice, the arrows go across the edges and the vertices 
are associated with the squares. There are 6 possible vertices, which can be divided into 
two types depicted in Fig.l, the vertices of each type being related by 7r/2-rotations. 
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(a) (b) 
Fig.l: The two types of vertices of the 6v model. 

For a special value of the spectral parameter the weights of the vertices of each class have 
the same Boltzmann weight. This specialization of the 6v model is known as the F-model. 
The weights of the two distinct vertices (a) and (b) of the F-model can be parametrized 
as 

(a) 0-6 = 1, (6) c = 2cos(7rA) (2.1) 
where the parameter A will be assumed to belong to the interval [0,1]. Following the 



standard argument of Baxter |jT^, we decompose the vertices into couples of line seg- 
ments making left or right turn as is shown in Fig. 2. The weights ( p.l| ) are then obtained 
by assigning a Boltzmann factor e^*'^/^ to each line segment, where the sign is positive 
(negative) if the line makes a left (right) turn. 
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(a) (b) 
Fig. 2: Decomposing the vertices as pairs of line segments. 

The sum over vertex configurations on each surface can be thus put in the form of a sum 
over fully packed closed nonintersecting loops. The Boltzmann weight of each loop C is 
the product of all phase factors along it. 



n 



2 ^ = e 



(2.2) 



where F is the geodesic curvature along the loop 

r{C) = ^ (# left turns - # right turns). (2.3) 

By the Gauss-Bonnet formula, on the flat lattice F = ±4 and the weight of each loop (after 
summing over both orientations) is 2cos(27rA). In the 6- vertex model one can introduce 
vortex operators of vorticity m, representing sources of m equally oriented lines. On the flat 
lattice m is restricted to be a multiple of 4. The 8-vertex model represents a deformation 
of the 6-vertex model obtained by adding vortices with m = ±4 with finite fugacity d. 

The 6-vertex model can be rewritten as a sohd-on-solid (SOS) model by assigning a 
variable ip G /32 to the sites of the dual lattice in such a way that the adjacent heights 
differ by ±/3/4. The oriented loops are then interpreted as domain walls of height ±/3/4. 
The m- vortex operator creates a discontinuity m/3/4 at a site. 

At large distances the height variable renormalizes to a free massless boson com- 



pactified on a circle of length /3 = 2nR |T^. If we choose the scale so that the duality 
transformation acts as 

1 47r2 

R^^ or P^- (2.4) 

then the conformal dimensions of the magnetic operators are 

^vortex _ ^vortex _ (w/3/27r) (2 5) 

Comparing the singularity of the free energy obtained from the exact solution of the 8- 
vertex model |lCm with the conformal dimension Al""^**^^, one finds 



(3 = TrVl -2A. (2.6) 

The Kosterlitz-Thouless point is at (3 = tt, where the lowest vortex operator (m = 4) 
becomes marginal. 

2.2. The 6-vertex model on a dynamical random lattice: loops, vortices, mapping to the 
compactified c = 1 string theory 

Since the weights of the F-model are invariant with respect to 7r/2-rotations, the 
model can be considered on any irregular square lattice S. The partition function can be 
again reformulated as a sum over loop configurations (Fig. 3), each loop C weighted by the 
geodesic curvature (|2.3| ) 

ns)= E n (2-7) 

loop conf loops C 
igurations 
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In this case the geodesic curvature (|2.3| ) depends on the geometry of the lattice through 
the gaussian curvature (the total deficit angle) TZ^C) enclosed by the loop C 

T{C) = 2n -n{C). (2.8) 

Therefore T{C) can be any integer (positive or negative) times 7r/2. 




Fig. 3: A piece of the random lattice covered by oriented loops. 



We define the partition function J-'(T, A, N) of the F model on a dynamical random 
lattice as the average of the partition function (|2.7|) in the ensemble of all connected 
abstract lattices obtained by gluing squares along their edges. The partition function 
depends on the "cosmological constant" T coupled to the area (^^ squares) A{S) and the 
"string coupling constant" 1/N coupled to the genus g{S) of the lattice: 



(2.9) 



It is natural to expect that the critical behavior of this statistical model is described by a 
compactified boson interacting with a Liouville field. Indeed, as it was shown in the 
exponential Z = e'^ can be identified with the partition function of a c = 1 string theory 
compactified at length = 27ri?, related to the vertex coupling as 



R 



1-A 



or 



/3 = 7r(l-A). 



(2.10) 



Here the scale is fixed again by adopting the convention ( p.4|) , according to which the 
self-dual radius is R = 1 {(3 = 2tv). For completeness, we give a short prove of the 
correspondence in Appendix A. Note that this is not the same correspondence as on the 
flat lattice and there is no reason that it should be the same. 

A vortex with vorticity m again represents a source of m equally oriented lines meeting 
at a point. However, since the lowest vortex charge is m = 1, an m- vortex operator creates 
a discontinuity m/3on the world sheet. It is necessarily accompanied by a conical singularity 
with curvature TZ = ^ (4 — m) . 

The scaling dimension of an m-vortex operator is obtained from the corresponding 
flat dimensionS ( p.5| ) by the KPZ rule |]15| and reads 



= ^ = (2.11) 

The scaling dimensions of vortex operators with m = 1, 2, 3, 4 are always smaller that one, 
which means that the interval < /3 < tt is deep in the vortex plasma phase. However, 
this is not a problem because in our model the fugacities of all vortices are strictly zero. 



3. The 6v matrix model 

3.1. From a complex matrix integral to a Coulomb gas system 



Ginsparg noticed in [0 that the exponential of the partition function (|2.9|) can be 
identified with the perturbative expansion of the following matrix integral 

f dXdxUxp (mrl-VrXX^ +X^X^^ - cos/3 (X^X)^]) 

^ e^(^'^'^) = , (3.1) 

j dXdX\e-^^^^^^^ 

where the integration goes over all x complex matrices X. 

At the point /3 = tt (A = 0), this matrix integral describes the dense phase of the 0{2) 
model on a random triangulation |]l|, or equivalently the compactified c = 1 string at the 
Kosterlitz-Thouless point /3 = tt. The correspondence with the compactified string for any 
13 has been worked out by S. Dalley [^. For reader's convenience we give the proof in 
Appendix A. 



^ Of course the length j3 is now given by ( 2.101) and not by ( ^.61) . 
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The vortex operators are introduced by adding a source term to the matrix potential, 
namely 

m 



^vortices = N ^{tvX^ + tr^t^ 



(3.2) 



m> 1 



The model ( |3.1| ) in presence of vertices (the random-surface analog of the 8- vertex model) 
has been solved in two particular cases: A = 0, di and d2 [|13| and c = 0, = b [|T6| . 

The complex matrix integral (|3.1| ) becomes gaussian at the price of introducing an 
auxiliary hermitian matrix A: 

/ X (3.3) 

W{X, A) = Ntr ( - + iAXX^e-'^/^ - iAX^Xe'^/^ - VTx^X j . 

The result of the integration over X is a matrix integral which can be viewed as a /3 7^ tt 
deformation of the 0(2) matrix model 



Z 



N 



z 



N 



dA 



det 



1 ® 1 - A®1 + ie'l^/^ 1 ® A) 



(3.4) 



The integrand depends only on the eigenvalues ai, ...,aN of the matrix A. After a liner 
change of variables, we write the partition function as 



Z 



N 



N 

n 

1=1 



-NV(ai) 



a,; 



n 



a,- — a 



where 



V{a) = - + 



a 



sm 



P 



In the following we will use the more convenient exponential parametrization 

a = -e"^ 
looks like 



in which the partition function 

r.00 N 



z 



N 



n 



-NV{<Pi) J-j- ^^^^ 2 



sinh i(^i — + i(3) 



(3.5) 



(3.6) 



(3.7) 



(3.^ 



where 



!>(</>) ^ V{-e'^). 



(3.9) 

Note that when is finite, the partition functions (|3.8| ) and ( |3.8| ) do not coincide, because 
the parametrization (|3.7|) assumes that a is negative. However, the two partition functions 
lead to the same espansion of the free energy and the observables. 
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3.2. The saddle point equations as difference equations for the resolvent 

We are interested in lattices with spherical topology, and therefore by the large 
limit of the integral (|3.8| ) . In this limit the spectral density 



1 ^ 

^(<^) = ]^E'^(<^-<^^) (3.10) 



i=l 



can be considered as a classical function, supported by some real interval [b, a] , where 
b < a. It is determined by the saddle point equation 



(3.11) 



tanh(^) tanh(^^^) tanh (^^^) 
For a time being we will consider a general potential of the form 

V(t)) = Y,tr.e^'' (3-12) 

n 

and later will specify 

T T 

h = -. (3.13) 



2sin^' ' 2 



Following the same logic as in [jT^, we will introduce the complex variat 

z^(p + it, (3.14) 
and formulate the saddle point equation (|3.11| ) as a set of conditions on the resolvent 

which is a holomorphic function on the cylinder z + 27ri = z cut along the interval [b, a] . 
The coefficients of the expansion of VF(z) at infinity 

W{z) = \ + Wie-' + W2e-^^ + ... (3.16) 



^ In the context of the 2D string theory, the complex z-plane represents the target space of the 
string. In terms of the c = 1 noncritical string, the real and imaginary parts of z are interpreted 
correspondingly as the Liouville and target-space directions. See, for example, the argument given 
in 
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are equal to the moments of the spectral density 

Wr.= r d(t>p{(l^)e^'^- (3.17) 
Jb 

The spectral density is equal to the discontinuity of the resolvent along the cut 

W{(j) + iO) -W{(j)-iO) = -27iip{(j)) (3.18) 
and the integral equation (|3.11| ) implies the following functional equation for the resolvent 



W{(f) + iO) + WU - iO) - WU + i(3) - WU - i(3) = ^^KM^ (3.19) 

dcp 

where it is assumed that (p G [6, a] . 

The normalization condition for the density can be written also as an integral along 
a contour C surrounding the interval [b, a] : 

1= rdM<P)= f^W{z). (3.20) 
Jb Jc 27rz 

3. 3. Geometrical meaning of the saddle point equations 
Let us introduce the new potential 



^W=E^^ (3-21) 



2 sinn-2 



related to the old one by 



and consider the meromorphic function 

j(z) = i[W{z + iP/2) - W{z - ip/2)] + U'{z). (3.23) 
The function J{z) is completely determined by the following four conditions: 

(1) - it is periodic with period 2z7r: J{z + ztt) = J{z — iir); 

(2) - it has two cuts {Rez G [6, a],Imz = ±i(3/2} in the strip \lmz\ < tt, and satisfies along 
the cuts a two-term functional equation 

J{(j) + i(3/2±iO) = J{<p-i(3/2TiO), (pe[b,a]; (3.24) 
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(3)- its series expansion at 2; — > +00 is 

j(z) = J^e^' + - 

^ ' sin/3 2sin/3/2 (3 25) 

+ 2sin(/5/2)W^i e"^ + 2 sin/3 VTae"^^ + ... 
where Wn are the moments of the spectral density defined by ( p.l7|) . 



(4)- it satisfies the normahzation condition ( |3.2C1| ): if Ci is a contour surrounding the upper 
cut, then 



/ ^Ji^) = ^- (3-26) 



The first two conditions mean that the function J{z) is a single valued analytic function 
on the torus obtained from the strip |Imz| < vr by cutting it along the intervals {Rez G 
[6, a],Im2 = ±i(3/2}, then identifying its two edges (Im^; = ±7r) as well as the "^^wei ^dge 
of the first cut with the ^^^^^r edge of the other one. The two main cycles of this torus are 
homotopic to the contour Ci surrounding the upper cut, and the contour C2 connecting 
the two cuts. 

The variable J is periodic along both cycles while the variable z is periodic along the 
cycle Ci and has discontinuity i(3 along the cycle €2- Therefore, if we find a canonical 
parametrization 

z = z{u), J = au) (3.27) 

for which the contours Ci and C2 are the images of the periods ui and UJ2 of a fiat rectangle, 
then the functions z{u) and C,{u) will satisfy the following periodicity conditions 

z{u + ui) = z{u), z{u + UJ2) = z + if], 

(3.28) 

C(w + Wl) = C(w), C{u + U2)=Ciu). 

Therefore, once the elliptic parametrization is found, the problem is essentially re- 
solvedi. The parameters of the solution are fixed by the normalization condition for the 
spectral density and the behavior of the function J{z) at z —>■ 00, which is determined by 
the potential U (z) . 

Let be the point in the periodic rectangle corresponding to z ^ 00. The map 
is assumed to be regular, which means that e~^*^'") can be Taylor expanded at u = Uoo- 
Therefore e^*^"-* ~ (m — Uoo)~^ and the third condition ( |3.25|) says that the function ((u) 
has a second-order pole at u = u^o and fixes the singular part of its Laurent expansion in 



Our method is a natural generalization of the method applied by J. Hoppe, V. Kazakov, the 



author and N. Nekrasov to resolve technically similar problems |19,2C], and which was originally 



proposed by J. Hoppe pl |. 
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3.4- Explicit solution in the elliptic parametrization 

We will identify the two periods of the torus with the standard quarter-periods of the 
Jacobi elliptic functions wi = 2K and U2 = 2iK' . The ^-cylinder with the two cuts is 
parametrized by the rectangle 

|Rew| < K, \lmu\ < K' . (3.29) 

The function C{z) is periodic with periods 2K and 2iK' and its only singularity in 
the rectangle (|3.29|) is a double pole aX u = Uoo- The residue of this pole is zero, since 
otherwise there would be another pole somewhere in the rectangle. The generic function 
with such properties (the Weierstrass elliptic function) depends on two constants A and B 



1 (3.30) 

an)=A + B^ -. 

sn^(u - Uoo) 

The function z{u) — is antisymmetric in u, which means that e^*-"-* has one simple pole 
at u = Uoo and one simple zero at u = —Uoo- It is therefore natural to look for it as a ratio 
of Jacobi theta functions. One can easily check that the function 

e^(-)= + (3.31) 

H{Uoo - u) 

satisfies the periodicity condition (|3.28| ) if the parameter Uoo is related to the period 2K 
by 

.oo = ^-^K = '-±^K. (3.32) 

The function ( p.31|) gives the desired regular map between the strip with two cuts in the 
z plane and the rectangle ( p.29| ), as is shown in Fig. 4. 

The two horizontal edges of the w-rectangle are mapped to the two cuts in the z-strip. 
Its two vertical edges are mapped to the line made of the segments [i^, in] and [— ztt, —i^]- 
The function z{u) along the "^owei given by 

z{±iK +y) = — — ± z- + In — -. (3.33) 

2 2 0{uoo-y) 

The endpoints of the cuts are the four branch points of the inverse map u = u{z), deter- 
mined by the condition dz/du = 0. The symmetry of the map ( p.31|) allows to determine 
the four branch points through a single parameter Ub G [0, K]: 

a±i(3/2 = z{ub±iK'), h ±i(3 /2 = z{-Ub ±iK'). (3.34) 
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The parameter Mf, is determined from 



c 



Ub) 



(3.35) 



where Z(u) 



eiu) 



Eiu) 



is the Jacobi Zeta function |j2^. In particular, if /3 = tt, 



then = Ub- FinaUy, the length of the cut is given, according to ( |3.34| ), by 

6(Woo + Ub) 



a — b — In 



OK 



Ub 



Note the following symmetry, of the solution, which reflects the symmetry (3 
(or A — > — A) of the original matrix integral 

(3 ^27T- 13, 
z{u) —z{K — w). 



(3.36) 

2n-p 




-K+iK' -"b+*' ' K+iK' 



-iK 



u 



K-iK' 



-K-iK' 



-iK' 



Fig. 4: The domains of the variables z and u. 
Now we have to fix, using the asymptotics (|3.25|) and the normalization condition ( |3.26|) , 
the four parameters of the solution, namely 



12 



- the coordinate of the middle point of the cut 

- the coefficients A and B, 

- the eUiptic no me /c^, which determines the two periods K and K' . 

The asymptotics at infinity gives three independent conditions, which together with 
the normahzation determine completely the four parameters. A relatively simple calcula- 
tion (see Appendix B) gives 

^_cot|^^ 
' " 2 H'{2u^y ^^-^^^ 



T* {K - E) [Z{2u^) + cndn(2t.^)] - K§j^ 



AT [zi2u^) + ^{2u 

Here we introduced the critical value of T 



2 

oo I 



(3.39) 



sin^^ 

T* = IGttA j- (3.40) 

cos ^ 

for which K' / K —>■ 0. The latter coincides with the critical coupling b* = 1/T* obtained 
by P. Zinn- Justin [Tl]. 



4. The scaling limit 

4.1. Explicit form of the solution in the scaling limit 

We would like to explore the scaling behavior near the critical point T ^ T* , where 

e^ 0, which in our parametrization corresponds to the limit K 00, K' ^ tv/2. In this 
limit the solution can be expanded in the dual modular parameter 

g = e-^/^'. (4.1) 

The expansion takes simpler form if we rescale slightly the variable u to 

u (4.2) 



2K 

which belongs to the domain |Imf | < ^, |Ref | < ^ In ^. 

13 



Our solution (|3.30|) - (|3.31| ) expands as a series in q as follows: 

1 ^--^ „ g^"" cosli2n(f — fc 



au)=A-B--B 



sinh^(f — f, 
sinli(t;oo + V 



+ 



n=l 



1-q 



2n 



=ln^^-^-(l + A)t; + ln . , , 
2 sinn(t; 



sinh 2nVoo sinh 2nv 



n 



1-q 



2n 



(4.3) 



(4.4) 



where 



1 + A, 1 
In -. 



2 g 

In order to find the q ^ asymptotics of ( [4.3|) and (f4.4| ) , we need the explicit expres- 
sions for the parameters T,a,b, B as functions of q. In the following we will use the symbol 
"pa" for "equal up to subleading terms in g". From ( |3.36| ), (|3.37|) , (|3.38| ) and (|3.39[) we get 



A 



T-T* 1 - A 1 



T* A ^ q 

a + b tan 1- 

g 2 ~ q * 



4A 



(1 + A)^(1-A) 

e 2 ^ ^ ^ - 



1- A 



T 

B ^ — 



2 

ttA 



T* (A + q1-^)2 
We see that indeed the position of the left cut tends to — oo 



(4.5) 



0, 



(4.6) 



while the right cut has a finite limit at g — > 



tan ^ , , ^ i+A , 

(1 + A)^(1-A) 



1- A 
2 



(4.7) 



Now we are ready to determine the scaling behavior of the solution (|4.3|) -( ^7^ ). The 
interesting scaling behavior is associated with the vicinity of the left pair of branch points 
of J{z), that is at 

b ± i— = z ( —Vb ± i—] ■ 
2 \ 2/ 

It is therefore convenient to introduce a new parameter r 

T = v + Vb, (4.8) 
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such that the two left branch points occur at r = ±|^. From ( p.35|) we get 

1, 1 + A 



oo -In^— . (4.9) 



and the point +oo of the 2-plane corresponds to 

Too ~ ^Voo - 2 I- \ ' 

Assuming that |r| ^ t^o and neglecting the subleading terms we get 



e^(") ^ C ((1 + A) e(i-^)" - (1 - A) e-(i+^)") (4.11) 

where _^ 

tan(7rA/2) 

(1 + A) 2 (1 - A) 4 

T/ie scaling behavior of the spectral density. Critical exponents associated with the 
entropy of the loops 

Let us first check that the "string susceptibility" exponent 7str is zero for all A. Con- 
sider the first momentum of the spectral density, whose singular behavior is that of a 
sphere with one puncture 

dlnZ 

Wi ~ ~ constant + (another constant) A^"''''^*^ 

The singular behavior is that of the coefficient S, which is given by the third equation 
( [1.5| ): B — B* ~ q^~^. On the other hand, by the first equation ( [1.5| ) 

gi-^ A 
A ~lni 

which implies that indeed 7str = and the free energy has logarithmic singularity 

InZ 



(4.13) 



Then we look for the scaling behavior of the spectral density and the resolvent. The 
spectral density is related to the function ({v) by 

p(.) = C(^ + W2)+C(«-W2) 
2n 
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and the resolvent ( p.l5|) is obtained by inverting the relation ( |3.23|) 

W{i7r + r) = W{i7r-T) = ^ C {^^ + i^) - C [-r - i^) ■ (4.15) 

It is easy to see that the function W{z) has only one cut on the physical sheet (the strip 
llm^l < tt), which goes along the interval z > b. On the other sheets there is an infinite 
number of cuts (at least for irrational A), placed at Imz = inn ± A, n G 2. 
At distances ^ e'', the relation between z and C, has two branches 

C±(^)~9^exp(^^^ (4.16) 

associated with the asymptotics at ^ 1. It is the smaller exponent which is the relevant 
one, but in the limit A ^ the two exponents coalesce and produce the observed logarith- 
mic behavior of the spectral density the 0{2) model. Correspondingly, the r-factor in the 
resolvent ( [1.1 5|) becomes in this limit(see eq. (3.41) of [^). The limit A ^ of the 



solution is worked out in Appendix C. 

Now let us find the scaling of an important observable: the typical length L of a loop 
on the world sheet. It is characterized by the critical exponent v 



L~v4^. (4.17) 

A related quantity is the average number M of loops on a worldsheet of area A. Since A 
is the total length of all loops, 

U=-^A^. (4.18) 
L 



To follow the analysis of [^J] and |2^, it is convenient to return to the original variable 
a = e'^ and consider the density 

p(a) = e-^p{4>) ~ aTTX (4.19) 
and the resolvent W{w) = (Tr ^^^ ) as a function of the complex variable w = —e^ 

W{w) = -e-^W{z)r^w^lnw. (4.20) 
The resolvent W{w) has a cut along the interval —oo<w<—M, where 

l + A 

A \ — 



A 



M = ~ ~ ( ) (4.21) 
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gives the scale for the length of the boundary. Since A ~ A and L ~ 1/M (for detailed 
arguments see [0), we conclude that 

ly = (4.22) 

1 + A ^ ' 

and the number of loops grows with the area of the world sheet as 

Af^A^. (4.23) 



Let us also mention that the exponent u determines also the fractal dimension of the e 
boundary of the disk amplitude, which is equal to 

Finally, let us calculate the dimensions of the m- vortex operators representing sources 
of m equally oriented lines. The correlation function of two such operators is obtained 
by calculating the "watermelon" configuration made of m nonintersecting lines relating 
the two points, as the one shown in Fig. 5. The singular behavior of such an amplitude is 
related to the scaling dimension S^'^^^^ of the m-vortex operator as 

Xm ~ A2'^'"-^^*^ (4.24) 




Fig. 5: A line configuration for the correlator of two 3- vortices. 



Proceeding in the same way as in Sect. 4 of ref.[2^|, one can extract the scaling di- 
mension of the m-vortex operator from the scaling behavior of the spectral density p(a) 
(eq. ( ^4.19| )). The amplitude corresponding to a configuration with m open lines is given 
(up to the phase factors e*^^ associated with the m lines, which can be neglected) by the 



convolutive integral in eq. (4.28) of ref. ||2J], and has the following singular behavior 



Compared it with ( [4.24] ), we find the same dimensions ( |2.11| ) 

1-A 



evertex 



m- 



(4.25) 



(4.26) 



which confirms the interpretation of the m-vortex operator as a discontinuity m/? on the 
world sheet. 
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5. Summary and discussion 

In this paper we gave the exact solution of the 6- vertex model on a dynamical random 
lattice with spherical topology. The model describes a gas of densely packed loops whose 
weights are coupled to the local curvature of the lattice. We have calculated the partition 
function of the disk with given length L and area A. For any value of the vertex coupling 
A, the critical singularity of the free energy corresponds to a theory with central charge 
c = 1, but the fractal dimension of the boundary of the disk changes continuously from 1/2 
to 1. The limiting value is achieved at A = 1, where the boundary has the same dimension 
as the bulk. This reminds what happens in the dense phase of the 0{n) model, where at 
n — s> the world sheet is all eaten by the boundary. Indeed, the exponent for the 
growth of the number of loops vanishes when A ^ 1, which means that in this limit the 
effective fugacity of the loops renormalizes to zero due to the fluctuations of the curvature. 
We also calculated the scaling dimensions of the m- vortex operators, representing sources 
of m equally oriented open lines on the world sheet. 

All our results agree with the interpretation of the 6-vertex model as a special real- 
ization of the compactified c = 1 string theory. From the string theory point of view the 
model is interesting with the possibility to treat very explicitly the vortex operators, which 
are difficult to construct in the standard realizations of the c = 1 string. 

The 6-vertex matrix model describes the c = 1 string theory compactified at length 
P < n, which is less than 1/4 of the Kosterlitz-Thouless length (3 = 47r, associated with the 
most relevant vortex operator (m = 1). In order to describe a theory compactified at length 
up to pn, it is sufficient to consider a chain of p coupled matrices. The diagonalization can 
be performed in a similar way and one obtains a chain of integral equations of the type 
( p.ll|) , which seems to be exactly solvable as well. 

The loop gas formulation of the c = 1 string theory places a bridge between the 



description with continuous target space via matrix quantum mechanics and the matrix 
models for the strings with discrete target space (the ADE models on a random lattice) 
||27|| , which describe all minimal conformal theories coupled to 2d gravity. In particular, 
it is clear why the string diagram technique based on a world sheet surgery works 
only for strings with discrete target space. A discrete target space without loops, i.e., a 
Dynkin graph of ADE type, can (and should) be immersed in the continuum by placing its 
adjacent points at distance /3 = tt, for which the curvature defects do not produce phase 
factors for the domain walls. In the case of the c = 1 string compactified at length P, this 
diagram technique works only at the point /3 = 1 (the 0{2) model). 
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The nonperturbative equivalence with the c = 1 string theory means that the 6v 
matrix model is exactly solvable in the so called double scaling limit . The integrability 
of this model follows also directly from the fact that the partition function (|3.8| ) is a r- 
function of the KP hierarchy. Using the scaling properties of the partition function integral, 
one can convert the KP equation into an ordinary differential equation, as it has been done 
inref. [0. 
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Appendix A. Equivalence to the c = 1 compactified string 

Let us integrate with respect to the angular variables in the original matrix integral 
( ^T|) . We decompose 

Xxt = UxU-^ 

with U G SU{N) and x = diag(xi, . . . , xat), and perform the U{N) integration with 
the help of the Harish Chandra- Itzykson- Zuber formula. Using the Cauchy identity, the 
integrand can be written as a determinant, and the partition function takes the form 



/■oo 

Zn = / dxi...dxN det K{xi,Xj) (A.l) 
Jo 



with 



K{x, x') = exp ( - ^nVt{x + x') + Nxx' - }rN{x^ + x'^) cos /?) . (A.2) 
We get rid of the linear term by a linear change in the integration variable 



y= v/iV^x-yo, yo= ^^^^'0^ =2^2NT/T\ (A.3) 

sm f 

after which the kernel becomes identical (up to normalization) to that of the inverse oscil- 
lator 

K{y,y') = —=== e ^.m^ ^ (^y\^P[,o^+ ,y )\y>^^ (A.4) 
A/iTT sm p 
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and the integration is to be done in the semi-infinite interval 



-yo <y < oc. 



Hence the partition function ( |A.1| ) describes a system of nonrelativistic fermions in 



inverse osciUator potential, stabilized by a wall at distance yo from the origina, at temper- 
ature /5 (see Fig. 6). The wall should be sufficiently far, so that the level of the Fermi sea 
is below the top of the potential. The critical distance is 



=2v/2(7r-/3)iV, (A.5) 

which corresponds to T = T* . Then the level of the Fermi sea reaches the top of the inverse 
gaussian potential. The string cosmological constant fi, which measures the deviation from 
the critical point, is proportional to {uq — yo)"^. 



U(y) 



i 













Fig. 6 : Fermi sea for the inverse gaussian potential 



The system of fermions on a circle of length P describes the singlet sector of the 
Matrix Quantum Mechanics (MQM) ||2^ and gives a non-perturbative realization of the 
c = 1 string theory compactified at length /3. The partition function of this system, 
calculated in ref. [^^, is invariant under the duality transformation 

P 271 

2^ ^ y ■ 

The equivalence with the compactified at euclidean time interval /3 is consistent only in 
the interval < /3 < Pkt, where Pkt = tt is the Kosterlitz-Thouless point. Near this 
point the the cutoff wall approaches the top of the inverse gaussian potential and there is 
no (metastable) ground state in the large limit. 



^ The restriction y > —yo is in fact imposed only on the trace, but it leads to a rapid exponential 
decay of the wave function of the fermions when y < —yo- 
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Appendix B. Fixing the parameters of the solution 

B.l. The three conditions following from the asymptotics at z ^ oo 

Using the explicit form ( |3.31D of z{u), we rewrite ( |3.25| ) as a series expansion in u—Uoo'- 



C{z) 



toe 



sin /3 



Hi2uc 



H'{0) J {u 



u 



g + b 



sin (3 



[H'iOW 



where 



tie^ H{2uoo) ^ e^+%H{2u^)H'{2u, 
2sinf H'{0) 

+ C 



T 



(B.l) 



2sinf 



fB.2) 



and the constant C is certain function of the elliptic modulus, whose explicit form will not 
be needed. This is to be compared with the first three terms in the Laurent expansion of 



B 



+ 



u 



^ 1 + o 

+ A-\ B 

3 

+ 0{Uoo - U). 



(B.3) 



which gives 



S-cot(/3/2) ^'(2^-^ 



8t2 H'^{2uc 



a + b 

e 2 



cotf H'{0) 



2 H'{2uoo) 
A^C-'-±^B. 



(B.4) 



B.2. The normalization condition 

The normalization (|3.20|) means that the integral of J{z) along the contour Ci sur- 
rounding its upper cut is equal to one. In terms of the variable u, this integral takes the 
form 

1 



27r 



K+iK' 



duz {u)C,{u) 



1. 



(B.5) 



K+iK' 
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Writing (using the standard notation Z{u) = 

z'{y + iK') = Z(uoo +y)+ Z{uoo - y) 

C{y + iK') =A + Bk%n^{uoo - y), 
we arrive at the following integral 



1 C 

1=7^/ dy[Z{uoo + y)+Z{uoo-y)][A + Bkhn\uoo-y)]. (B.6) 



With the help of the addition formula [31], 17.4.35, we write it in the form 



1 = -—[Z{2u^)h + sn{2u^)l2] 

where „ 

/•^ ^ 2 2{K-E) 
h = av sii V = 

J-K ^ 

I2 = / dv sn(2uoo — v) sn^v. 

J-K 

The second integral can be decomposed as 

cndn 
sn 

where 



K 2 

, sn^t; 
dv- 



_x 1 — /c^sn2(2Moo)sn2t) 
is a standard elliptic integral 435.02 : 

2K Z{2u^) 



I 



2 



/c2sncndn(2uoo) 
Substituting this in (|B.6| ) yields 



^ = 2(K - E)[Z{2u^) + cndn(2^.^)] - 2^^^. (B.7) 
B sn(2uoo) 

Appendix C. The special case (3 = n. 

Here we check that the solution of the 0{2) matrix model with quadratic potential is 
obtained as the limit /3 ^ tt (or A ^ 0) of our general solution ( |3.3U| )-( |01] ). This limit is 
quite subtle and does not commute with the scaling limit q ^ Q. In this limit we have 



T* = 32, — = ^4-, e« = — , = k'e\ (C.l) 

^ 8E ^ ^ 
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At A = we have more symmetry because the special points of the map coincide with the 
quarter of the period, Uoo = ut = K/2, and the function e^*^"-* becomes a Jacobi eUiptic 
function 



e^(") = dn - y - iK' 

The function ({u) is analytic in 2; at A = 

B 









/ sn 





(C.2) 



C(^^) =A + 



but the coefficient B tends to infinity as 1/A 

TtT/T 



A + B - Be-^'^+^^e^^^''^ 



(C.3) 



B 



A ■ 



(C.4) 



Therefore, in order to take the limit A ^ 0, we should first subtract the diverging analytic 
piece of (. The finite, nonanalytic piece is then given by 



Ciu) = BX 



d 



where 

Further we write 
d 



dX sn2(w - K/2 - XK/2) J 



d _ d f du\ d 
dX^dX^ \dx)^d^' 



(C.5) 



1 



A=0 



dX sn^{u- K/2 - XK/2) 

K-2{dxu)^ = 2K 



{K - 2{dy,U), 



cndn 



sn" 



K 



u — 



H' sncn 



H dn 

Since k' [|^] (w + -^) = — e^~", we can write C,{u) 

KT 



K 



as 



an)- 

The corresponding spectral density 



2z(u-) d , K\ 
-e-(^)-lniy^. + -j. 



(C.6) 



C{u + iK') + C{u - iK') 
2^r 



coincides with the one obtained by solving directly the saddle point equations for the 0(2) 
model [02 . 
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